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In this work we present results for pairing gaps in /3-stable neutron star matter with electrons and 
muons using a relativistic Dirac-Brueckner-Hartree-Fock approach, starting with modern meson- 
exchange models for the nucleon-nucleon interaction. Results are given for superconducting ^So 
protons and ^P2 and neutron superfluids. A comparison is made with recent non-relativistic 
calculations and the implications for neutron star cooling are discussed. 

PACS number(s): 97.60. Jd 21.65.+f 74.25.Bt 

Superfluidity and superconductivity of matter in neutron stars is expected to have a number of consequences directly 
related to observation, see Refs. Among processes that will be affected are the emission of neutrinos. Neutrino 

emission from e.g. various URCA processes are expected to be the dominant cooling mechanism in neutron stars less 
than 10^ — 10^ years old. Typically, proton superconductivity reduces considerably the energy losses in so-called 
modified URCA processes and may have important consequences for the cooling of young neutron stars. Another 
possible manifestation of superfluid phenomena in neutron stars is glitches in rotational rates observed in a number of 
pulsars. Moreover, the estimation of superfluid gaps and studies of pairing are not only important issues in neutron 
star matter, but also in the rapidly developingfield of neutron-rich systems such as heavy nuclei close to the neutron 
drip line |^ or the study of light halo nuclei [0] . Therefore, theoretical studies of pairing in neutron-rich assemblies 
form currently a central issue in nuclear physics and nuclear astrophysics. 

The aim of this Letter is to present results from self-consistent calculations for neutron and proton pairing gaps in 
/?-stable matter relevant for neutron star studies. Pairing in the partial waves ^5*0, "^^2 a-nd will be studied 
using a relativistic approach with modern meson-exchange potential models to describe the nucleon-nucleon (NN) 
potential. A comparison with the corresponding non-relativistic approach is also made. 
Our computational scheme is as follows: 

The first ingredient in our calculation is the self-consistent evaluation of single-particle energies in /3-stable matter 
starting from the meson-exchange potential models of the Bonn group 1^ . These single-particle energies are obtained 
within the framework of the Dirac-Brueckner-Hartree-Fock (DBHF) scheme [PHl l| , using a medium renormalized 
NN potential G defined through the solution of the G-matrix equation 

GM = y + yQ — l^QGH, (i) 

uj - QHoQ 

where u is the unperturbed energy of the interacting nucleons, V is the free NN potential, Hq is the unperturbed 
energy of the intermediate scattering states, and Q is the Pauli operator which prevents scattering into occupied 
states. Only ladder diagrams with two-particle states are included in Eq. (|^). In this work we solve Eq. using 



the Bonn A potential defined in Table A. 2 of Ref. This potential model employs the Thompson |10 IJ] reduction 
of the Bethe-Salpeter equation, and is tailored for relativistic nuclear structure calculations. For further details, see 
Refs. 

The DBHF is a variational procedure where the single-particle energies are obtained through an iterative self- 
consistency scheme. To obtain the relativistic single-particle energies we solve the Dirac equation for a nucleon in the 
nuclear medium, with c = h ^ 1 , 

{^-m + ^ip))uip,s)^0, (2) 

where m is the free nucleon mass and u(p, s) is the Dirac spinor for positive energy solutions, with p — {p'^,p) being 
a four momentum and s the spin projection. The self-energy for nucleons can be written as 

- Ssb) - 7oS"(p) + 7PS^(p). (3) 

Since << 1 we approximate the self-energy by 

S « Ss - 7oS° = Us + Uv, (4) 

where Us is an attractive scalar field and Uy is the time-like component of a repulsive vector field. The Dirac spinor 
reads then 
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where Xs is the PauU spinor and terms with tilde Uke Ep — -^p^ + represent medium modified quantities. Here 
we have defined |l0| , [l^ rh — m + Us- In all equations below, a momentum p refers to the three-momentum p. The 
single-particle energies ip can then be written as 

ep = Ep + Uv (6) 
The single-particle potential Up is given by Up — Usifi/Ep + Uv and can in turn be defined in terms of the G-matrix 

eTTT ^P^\ '^(^ ^£p + ^h) \ph) , (7) 

h<kj, ^hEp 

where ph represent quantum numbers like momentum, spin, isospin projection etc of the different single-particle states 
and kp is the Fermi momentum. Eqs. are solved self-consistently starting with adequate values for the scalar 

and vector components Us and C/y The proton fraction in /3-equilibrium is determined by imposing the relevant 
equilibrium conditions on the processes e~ + p — > n + i^e and e~ — > jjT +v^ + Ve. The conditions for /3-cquilibrium 
require that /i„ = fip + fie, where Hi is the chemical potential of particle species i , and that charge is conserved 
^p = rie, where is the particle number density for particle species i. We also include muons and the condition for 
charge conservation becomes Up = rie + n^, and chemical equilibrium gives fj,e — fJ-^,- Throughout we have assumed 
that neutrinos escape freely from the neutron star. The proton and neutron chemical potentials are determined from 
the energy per baryon, calculated self-consistently in the above DBHF approach. 

The next step in our calculations is to evaluate the pairing gaps for various partial waves. To evaluate the pairing gap 
we follow the scheme of Baldo et al. jlj] , originally proposed by Anderson and Morel . These authors introduced 
an effective interaction Vk^k' ■ This effective interaction sums up all two-particle excitations above a cutoff momentum 
fcjvf, = 3 fm""'^ in this work. It is defined according to 



Vk^k' =Vk,k' - ^ Vk^k" -^^Vk" ,k' , (8) 

k">k 



where the energy £k is given by £k = y {ik ^ £f) + ep being the single-particle energy at the Fermi surface, 
Vjt,fe' is the free nucleon-nucleon potential in momentum space, defined by the three-momenta fc, k' . The renormalized 
potential Vk,k' and the free NN potential Vk,k' carry a factor / EkEk' , due to the normalization chosen for the Dirac 
spinors in nuclear matter. These constants are also included in the evaluation of the G-matrix, as discussed in pO| , pl| . 
For the ^5*0 channel, the pairing gap is [|6j 

A. = - ^ Vk,k'^. (9) 

For the partial wave we employ the expressions given in Rcf. ||l7[| , modified as well by the above normalization 
constants. For further details, see e.g. Refs. In summary, first we obtain the self-consistent DBHF single- 

particle spectrum for protons and neutrons in /3-stable matter using the method detailed in Ref. fl^ . Thereafter 
we solve self-consistently Eqs. and in order to obtain the pairing gap A for protons and neutrons for different 
partial waves. 

Our results for the pairing gaps, scalar and vector potentials for neutrons and protons, proton and neutron fractions 
and the chemical potential for electrons (and muons for total baryonic densities greater than p — 0.15 fm^'^) are 
displayed in Tables | and || as functions of the total baryonic density. The results of these tables can in turn be 
used in relativistic equations for various modified URCA processes, in a similar way as done in the non-relativistic 
approach of Friman and Maxwell ||20|] . In Fig. |l] we plot as function of the total baryonic density the pairing gap 
for protons in the ^Sq state, together with the results from the non-relativistic approach discussed in Refs. p8| , p^ . 
The results in the latter references were also obtained with the Bonn A potential of Ref. These results are all 
for matter in /3-equilibrium. In Fig. |^ we plot the corresponding relativistic results for neutron energy gaps in the 
■^^2 channel. For the ^Z?2 channel we found both the non-relativistic and the relativistic energy gaps to vanish. The 
non-relativistic results for the Bonn A potential are taken from Rcf. [^. We have omitted a discussion on neutron 
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pairing gaps in the ^Sq channel, since these appear at densities corresponding to the crust of the neutron star. The 
gap in the crustal material is unlikely to have any significant effect on cooling processes 

As can be seen from Fig. ^ there are only small differences (except for higher densities) between the non-relativistic 
and relativistic proton gaps in the -'^iS'o wave. This is expected, since the proton fractions (and their respective Fermi 
momenta) are rather small, see Table |[ For neutrons however, see Table ||, the Fer mi momenta are larger, and we 
would expect relativistic effects to be important. At Fermi momenta which correspond to the saturation point of 
nuclear matter, kp — 1-36 fm~^ the lowest relativistic correction to the kinetic energy per particle is of the order 
of 2 MeV. At densities higher than the saturation point, relativistic effects should be even more important, as can 
clearly be seen in the calculations of Ref . . Since we are dealing with very small proton fractions in Table ||, a 
Fermi momentum of kp = 1-36 fm~^, would correspond to a total baryonic density ~ 0.09 fm"'^. Thus, at larger 
densities relativistic effects for neutrons should be important. This is also reflected in Fig. |^ for the pairing gap in 
the ■^P2 channel The relativistic gap is more than half less than the corresponding non-relativistic one, and the 
density region is also much smaller. This is mainly due to the inclusion of relativistic single-particle energies in the 
energy denominator of Eq. ^ and the normalization factors for the Dirac spinors in the NN potential. Even the 
non-relativistic energy gaps for neutron star matter in /3-equilibrium are small compared with the results for pure 
neutron matter, where the '^^2 energy gap has a maximum around 0.12 — 0.13 MeV, see Refs. [0,|l]j22). The 



consequences for cooling rates and the interior composition of a neutron star are significant. A recent investigation 
of various cooling mechanisms by Schaab et al. p3t found that an agreement with observed surface temperatures was 
reached if the ■^P2 energy gaps were of the order 0.05 MeV. Our non-relativistic results for /3-stable matter are of 
this size, while the relativistic energy gaps result in an almost negligible suppression of e.g. various modified URCA 
processes in the interior of a neutron star. These results, and those of Schaab et al. as well, are at askance with 
those of Page Q|, where, in order to explain the observed temperature of Geminga, baryon pairing has to be present 
in most, if not all of the core of the star. 

In summary, in this work we have calculated in a self-consistent way single-particle energies and energy gaps using a 
relativistic DBHF approach. To our knowledge, after the relativistic work of Kucharek and Ring ||l6| , this is the first 
estimate of pairing gaps within the framework of the DBHF approach. In Ref. ||l6|l , the -"^^o gap in symmetric nuclear 
matter was studied within the framework of the Serot-Walecka model [ p^ . The only parameters which enter our 
approach are those which define the free NN potential Here we have focused on pairing in dense matter, though 
our approach allows also for a consistent treatment of other neutron star properties. The same NN force used here 
has also been used in Ref. |2j] to calculate the equation of state and the total mass and radius for a neutron star. 
Combining the results from this work and those of Refs. ||l^,Q, the following picture emerges: 

Within the DBHF approach, the direct URCA processes are only allowed for densities larger than 0.52 fm""^, see 
Ref. ijl^. A neutron star with total mass I.6M0 would have a central density of pc — 0.4 fm^'^ in /3-stable matter 
For such a central density, various modified URCA processes are possible mechanisms for neutrino production 
in a neutron star. The main suppression of these processes would then come from protons in the ^Sq state. The 
reader should note that there are other possible cooling mechanisms than those discussed here, such as neutrino-pair 
bremstrahlung p5|] , direct URCA with hyperons or neutrino emissions from more exotic states, such as pion and kaon 
condesates or quark matter, see e.g. Refs. [^|j5p^ for recent reviews. However, for a star with central density pc = 0.4 
fm""^, many of these more exotic neutrino emissivities are less likely. Hyperons appear at densities p ^ 0.3 or greater 
p6t . Similar densities are expected for kaons and quark matter ]23| , p6t . In addition, neutrino-pair bremstrahlung was 
recently found |25| to be much less important than previously estimated. Thus, for a 1.6Mq neutron star with central 
density of 0.4 fm^"^ obtained with our DBHF approach |24j, the most likely cooling scenario is through modified 
URCA processes, and the main suppression comes from superconducting protons in the ^S'o state. Finally, it ought 
to be noted that we have not included effects from medium polarization effects, as discussed in Ref. These may 
further change the size of the energy gaps and the neutrino emissivities (the ^5*0 gap should decrease while the 
gap is expected to increase). Further, we have not considered the possibility of pairing, which appear due to the 
increased proton fraction, as discussed by Aim et al. p^ . 

This work has received support from The Research Council of Norway (NFR) (Programme for Supercomputing) 
through a grant of computing time. MHJ thanks the Istituto Trentino di Cultura, Italy, and the NFR for financial 
support. 
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TABLE L Proton fractions Xp, scalar and vector single-particle potentials and Uy, respectively, for protons, the proton 
pairing gap Ap for protons in the ^ So state and the electron (and muon) chemical potential as functions of total baryonic 
density p. Densities are in units of fm~'^ , [/|, Uy, Ap and /ie in units of MeV. 



p 


Xp 


ui 




Ap 


Me 


.0013 


.0032 


-7.8479 


3.2471 


.0121 


11.7231 


.0068 


.0050 


-77.7002 


61.7252 


.0483 


20.3904 


.0281 


.0096 


-172.0541 


135.3744 


.2024 


38.9884 


.0583 


.0156 


-236.5725 


181.5207 


.4386 


58.5459 


.0944 


.0229 


-285.0128 


213.1141 


.7036 


78.1881 


.1377 


.0307 


-329.1642 


242.7944 


.9107 


98.3550 


.1811 


.0403 


-365.8355 


270.4411 


1.0160 


115.8907 


.2007 


.0462 


-381.3338 


283.5829 


1.0173 


123.0215 


.2212 


.0524 


-396.7707 


297.3635 


.9742 


130.1985 


.2627 


.0658 


-424.5634 


325.2710 


.7712 


143.9456 


.3072 


.0801 


-451.9637 


357.1098 


.4490 


158.2441 


.3304 


.0877 


-464.7640 


373.9551 


.2638 


165.5386 


.3544 


.0953 


-476.8407 


391.2967 


.1826 


172.9228 


.3594 


.0968 


-479.2122 


394.8924 


.0856 


174.4599 
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TABLE II. Proton fraction Xp, neutron scalar and vector single-particle potentials Ug and Uy, respectively and the neutron 
pairing gap A(^P2) as functions of total baryonic density p. Densities are in units of fm~^ , Ug, U" and A in units of MeV. 









'■v 




.0756 


.0191 


-118.4076 


90.1259 


0.009 


.0811 


.0202 


-127.8562 


97.9057 


0.013 


.0849 


.0210 


-134.2159 


103.1913 


0.014 


.0949 


.0230 


-150.7538 


116.9925 


0.017 


.1012 


.0243 


-161.1272 


125.6867 


0.017 


.1056 


.0252 


-167.9521 


131.2468 


0.017 


.1125 


.0266 


-179.0345 


140.6626 


0.015 


.1172 


.0275 


-186.6448 


147.1867 


0.013 


.1196 


.0279 


-190.5106 


150.5173 


0.011 
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o : relativistic 

X : non-relativistic 
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FIG. 1. Proton pairing in /3-stable matter for the ^5*0 partial wave. The non-relativistic results are taken from Ref. [18]. 
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FIG. 2. Neutron pairing in /3-stable matter for the P2 partial wave. The non-relativistic results are taken from Ref. [17]. 
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